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We irradiate the heavy fermion superconductors (U,Th)Bei3 with high-energy heavy ions. Damage from the 
ions affects both heat capacity and magnetization measurements, although much less dramatically than in other 
superconductors. From these data and from direct imaging, we conclude that the irradiation does not create the 
amorphous columnar defects observed in high-temperature superconductors and other materials. We also find 
that the damage suppresses the two superconducting transitions of Uo.grTho.osBeia by comparable amounts, 
unlike their response to other types of defects. 
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I. INTRODUCTION 

Microscopic lattice disorder can have interesting effects, 
particularly on unconventional superconductors. Sensitiv- 
ity of the critical temperature Tc to non-magnetic scatterers 
serves as evidence of non-s-wave electron-pairing. Defects 
can be introduced during sample fabrication, or induced af- 
terwards by irradiation. The latter possibility encompasses 
both point defects created by light ions and other particles, 
and amorphous columns created by high-energy heavy ions. 
These columnar defects, which have been heavily studied in 
high-temperature superconductors (HTS), serve as highly effi- 
cient pinning centers for vortices and alter a sample's response 
to a magnetic field. 

Little work has been done on another group of unconven- 
tional superconductors, the heavy fermion (HF) systems. The 
two orders of magnitude difference in the critical temperatures 
of HTS and HF superconductors changes the scale for thermal 
effects. Thermal fluctuations which complicate HTS behavior 
should be absent in HF. 

The spectacular response of one heavy fermion material, 
UBei3, to thorium doping makes it a good candidate for in- 
vestigating other forms of lattice disorder. Tc{x) is not mono- 
tonic for Ui_2:Th2;Bei3, and a second thermodynamic phase 
transition appears for 0.019 < x < 0.043 |1, 2]. No other 
dopant has such effects. The Be as well as the U is sensi- 
tive to substitution details. The superconducting transition is 
suppressed below 0.015 K for UBe12.94Cuo.06 1 3], while the 
same amount of B has virtually no effect |4]. 

Multiple superconducting phases indicate an unconven- 
tional order parameter. A second HF material, UPts, also ex- 
hibits a double transition |I5]. Furthermore, the lower phases 
of both UPt3 and Ui_:EThj;Bei3 have a significantly reduced 
vortex relaxation rate, while in pure UBei3 the relaxation rate 
is linear in temperature 0]. This may be evidence of time re- 
versal symmetry breaking in the lower phase. If time reversal 
symmetry is broken, spontaneous magnetic fields can appear. 
Boundaries between different magnetic domains can obstruct 
vortex motion, leading to the observed drop in relaxation rate. 



Vortex motion is also reduced in Ui_a:Th2;Bei3 for zero-field- 
cooling as compared to the field-cooling case |7]. Domain 
wall pinning may explain this effect as well, since cooling in 
a field should lead to larger domains and fewer domain walls. 
The atypical vortex behavior and the distinction between the 
two superconducting phases suggests that vortices may also 
have interesting interactions with other types of pin sites, such 
as columnar defects. 

Here we investigate heavy-ion irradiation of (U,Th)Bei3. 
We show that calculations predict columnar tracks to form. 
However, heat capacity measurements, imaging, and magneti- 
zation measurements suggest that such tracks are not actually 
present. We also discuss how heavy-ion damage affects the 
different superconducting phases. 



II. IRRADIATION DAMAGE 

The thermal spike model explains damage creation in both 
elements and alloys |8, 9, 10]. In this model, the heavy ion 
loses energy primarily to the electrons, which in turn ex- 
cite phonons. The electron-phonon coupling strength g de- 
termines how fast the energy is transfered to the lattice, while 
the thermal conductivities and Kp of electrons and phonons 
govern how the energy spreads through the sample. The cou- 
pled differential equations 



c — - —(k —) 

^ dt dr dr r dr 



g{T,-Tp) + A{r,t) 



d dTp Kp dTp 



dr dr 



dr 



g{Te - Tp 



describe the energy transfer. Here Kg and the specific heat 
Ce of the electron system are functions of the electron tem- 
perature Te, while Kp and the lattice specific heat Cp depend 
on the phonon temperature Tp. Se = ^^|e is the energy 
loss from the heavy ion to the electrons in the target. As Se 
increases, a threshold value allows formation of amorphous 
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TABLE I: Parameters for thermal spike calculations in UBeis. 
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tracks. With further increase the track radius grows. For the 
radial distribution of the original ion energy loss we use 
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asinCUl- Herein = 19/31/3, L = 0.1 nm, M = (1.5 + 0.5^) 
nm, and (3 is the ion velocity relative to the speed of light. The 
maximum range r,,nax = 6 x 10~^( ^^^^ 1.079 ^j^^^ where m 
is the electron mass in grams and c the speed of light in cm/s. 
The time to is of order lO^^^ s and 6* = 9.84 x 10"^ cm. 
The normalization of A{r, t) is chosen so that the total energy 
transfer to the electrons equals Se- This formula, particularly 
the term involving K, L, and M, comes from fits to observed 
radiation distributions in water 1 11]. Thermal spike calcula- 
tions using these equations successfully predict the existence 
and even radius of columnar defects in a variety of pure metals 
and alloys 1 8 ] . 

Numerical solution of these equations requires knowledge 
of various parameters of the material. Table U shows the val- 
ues we use. Since most of the simulation is at high temper- 
ature, we use high- temperature limits in several cases. The 
electronic and lattice specific heats are Cg = l.SfcsTig and 
Cp — iksna- The density of atoms in the lattice is Ua, while 
the electron density is = zua with z the average valence 
of the atoms composing the material. For the electronic ther- 
mal conductivity Kg we use the Wiedemann-Franz law. Al- 
though in most metals the high-temperature electrical resis- 
tivity is proportional to temperature, in UBei3 the resistivity 
is nearly constant, so Ke is itself proportional to temperature 
fT?]. For the lattice thermal conductivity Kp we expect a 1/T 
temperature dependence and obtain the prefactor from the 90 
K value O. 

The heat of fusion of a compound equals the latent heat of 
its component elements, plus the difference between the liq- 
uid state mixing enthalpy and the solid enthalpy of formation. 
Typical values are 10^° to 10^^ erg/cm"^. In this case the U 
and Be latent heats contribute 2 x 10^" erg/cm^^. The heat of 
formation at room temperature is —2 x 10^° erg/cm'' 1 13], so 
4 x 10^° erg/cm'3 is an upper bound on the heat of fusion. Al- 




FIG. 1: Predicted defect radius in UBeis. The two curves use heat of 
fusion 2 X 10^° erg/cm^ (•) and 4 x lO"' erg/cm^ 
simulation of electronic energy loss for 238u67+ 
(•) ions in UBei3. The arrows indicate the energy of the beams used 
for our irradiations. 
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though the mixing enthalpy is unknown, 2 x 10^*^ erg/cm^ is 
a more reasonable estimate for the heat of fusion. 
Finally, for the electron-phonon coupling we use 0] 
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Here L — 0.245 erg-fi/K is the Lorentz number 

Figure ^ shows the defect radius obtained by solving the 
thermal spike equations numerically. The values are compara- 
ble to those of Ce and Y alloys ] 9], which have similar resistiv- 
ity, Debye temperature, and heat of formation to UBeia. The 
high atomic density of UBei3 increases its sensitivity to ir- 
radiation, while the unusual temperature independence of the 
resistivity decreases the sensitivity. 



III. SAMPLE PREPARATION AND IMAGING 

Our samples are polycrystalline UBei3 and 
Uo.97Tho.o3Bei3 pieces with a grain size of about 100 
/im, prepared by arc melting in an argon atmosphere and 
subsequent annealing at 1400° C in a beryllium atmosphere 
for 1000 hours ] 17]. Before irradiation, samples are thinned 
to 25 /im to prevent ion implantation, with a typical lateral 
dimension of 1 mm. The irradiation is performed at ATLAS, 
Argonne National Laboratory with 1.4 GeV 238^67+ 
208PI356+ j^jjg jj^g Jqjj velocities are 5.9 MeV/amu and 6.7 
MeV/amu, respectively. The matching fields _B$, at which 
the vortex density and track density are equal, range from 0. 1 
T and 7 T, corresponding to typical lateral track separations 
from 155 to 18.5 nm. The 208p^56+ -^^^^ ^^^.^ ^^^^ ^.j^^ 
3 T and 7 T samples, while the 238^67+ ^gj.g \i^,e6. for the 
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TABLE II: Effects of irradiation on transitions in (U,Th)Bei3. 
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Transition temperatures marked by (*) come from magnetization measurements, the remainder from heat capacity. 



remaining samples. The irradiation takes up to two hours, 
depending on dosage. Because of the very small probabiUty 
that two ions will reach the sample at nearly the same time 
and position, we assume that each ion acts independently in 
damaging the sample. 

The inset of Figure ^ shows the electronic energy loss Se 
as a function of the incident beam energy, obtained from a 
Stopping and Range of Ions in Matter (SRIM-98 code lUsjl l 
simulation. The beam energy used for the irradiation is near 
the maximum and gives Se = 31.2 keV/nm for uranium, 
Se = 27.6 keV/nm for lead. The corresponding nuclear en- 
ergy losses are negligible, 0.05 keV/nm and 0.04 keV/nm, re- 
spectively. The predicted stopping distances are 56 /im and 60 
/im, more than twice the sample thickness. From our thermal 
spike calculation, the irradiation should produce defects with 
radius 2.8 nm (uranium) or 2.2 nm (lead). 

After irradiation, both the 2.0 T UBeia and 7.0 T 
Uo.gvTho.osBeia crystals were further thinned and viewed 
in both a Philips CM-30 transmission electron microscope 
(TEM) and a JEOL 3000F high resolution transmission elec- 
tron microscope. The planes used for the images on the JEOL 
machine limited the resolution to about 0.5 nm. At _B$ =2.0 
T the defects should be about 32 nm apart. The total viewing 
area of about 73,000 nm^ should contain about 70 defects but 
in fact shows no evidence of amorphous columns. For the 7.0 
T sample, the imaging area of over 54,000 nm^ should have 
nearly 200 defects. We again find no damage of radius greater 
than 2 nm, and only six candidates with radius from 1 nm to 
2 nm. Of these, closer examination shows that three are edge 
dislocations. The remaining candidates could also be dislo- 
cations, viewed from angles that do not clearly show a lattice 
plane ending. We find no compelling reason to identify them 
with columnar tracks, since they do not match the expected 
defects either in size or in density. We conclude that colum- 
nar defects are either absent in our sample or far smaller than 
anticipated from our thermal spike model calculations. 

We emphasize that the microscopes we used can detect 
amorphous regions of the expected size. A transmission elec- 
tron microscope successfully imaged columnar defects of ra- 
dius 5 nm in UO2 [19]. High resolution TEM has detected 
columnar defects with radius down to J^nm, for example in 
tin oxide |j20| and BiaSraCaCuzOs+a ||2ljl. 



IV. HEAT CAPACITY MEASUREMENTS 

We measure heat capacity on a dilution refrigerator at tem- 
peratures down to 100 mK and magnetic fields up to 8 Tesla. 
We use a relaxation method with a RUO2 thin film thermome- 
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FIG. 2: Zero-field heat capacity of unirradiated (□), 5 T irradiated 
(•), and 7 T irradiated (A) Uo.grTho.osBeis. 



ter, a [50:50] AuCr thin fikn heater, and a fine copper wire as 
a heat link. 

Figure 12 shows C/T for the unirradiated, 5 T, and 7 T 
Uo.gvTho.osBeia samples. Since we do not know the precise 
sample sizes, we use a multiplicative constant to set the nor- 
mal state value of C/T to 1.1 J/mole for each curve. We 
do not adjust the measured heat capacity for any contribution 
from the thermometer, heater, or mount; but comparing the 
curves shown to previous heat capacity measurements on bulk 
samples suggests that background effects are less than 30% of 
our signal. In any case, the background should not affect our 
determination of the Te suppression or transition width. 

Irradiation does not measurably change Tc up to 5$ = 2 
T. As shown in Figure|2]for Uo.gyTho. 036613, C/T changes 
shape only slightly even at a defect density of 5 T, broadening 
by several mK [ 22 1 . At the higher defect density of = 7 T, 
the transitions do change shape. The upper transition widens 
by roughly a factor of two, and both decrease significantly in 
amplitude. Here width is defined by identifying the change 
in C/T across the transition and using the temperature dif- 
ference between the 10% and 90% values of the heat capacity. 
Irradiation has a similar effect on pure UBeia, reducing Tc and 
slightly broadening the transition, while retaining the shape of 
the C/T curve up to 5 T. TableHIIsumm arizes the effects of ir- 
radiation on transition temperatures for both compounds. The 
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marked transitions come from magnetization measurements 
rather than heat capacity, as will be discussed later 

Since the 208p|^56+ j^j^^ used for irradiation have less en- 
ergy loss than the 238^67+ j^jj^^ ^jjgy produce smaller defects 
in a thermal spike calculation. Depending on the value of the 
heat of formation, the uranium atoms might even produce de- 
fects while the lead atoms do not. We note here that the Tc 
depression for the 3 T sample is consistent with the others. 
The 7 T samples have, if anything, a larger change in than 
expected from the lower irradiation dosages. Thus it appears 
that the 208p^56+ ^^^^ ^^^^j^ damage as the 238^67+ 

ions, despite their smaller energy loss. 

It is tempting to take the nearly identical percentage sup- 
pression of the two thoriated transitions as evidence that 
both superconducting phases have the same pairing symmetry. 
However, the lower transition is between two superconduct- 
ing phases, rather than between a superconducting and a nor- 
mal phase. The usual arguments on how impurities affect the 
transition do not apply. We do note that the heavy-ion irradi- 
ation produces a different response from either non-magnetic 
or magnetic substitutional impurities |23]. Non-magnetic (La) 
impurities depress the upper transition but leave the lower 
transition unchanged, while magnetic (Gd) impurities sup- 
press the lower transition more than twice as strongly as the 
upper one |23]. These contrasting behaviors of the two tran- 
sitions for different types of defects emphasize the extreme 
sensitivity of the (U,Th)Bei3 system to changes in electronic 
structure. 

The Tc suppression proves that the irradiation did create 
some damage within the crystals. However, the small amount 
of suppression, like the lack of visible columns in the imag- 
ing, suggests that the damage is not in the form of amorphous 
columns. For comparison, U ions at a 2 T matching field re- 
duce Tc of YBa2Cu307_a; by 5.3% 124], more than we find at 
5 T. Although in principle our low Tc reduction could occur 
from insensitivity to defects, rather than an absence of defects, 
HF superconductors are more sensitive than other supercon- 
ductors to both neutron and light-ion irradiation [ 3, 25]. Even 
our Bcf, — 7 Tesla irradiation affects Tc far less than the point 
defect work 



V. MAGNETIZATION MEASUREMENTS 

The motion of superconducting vortices is particularly sen- 
sitive to defects. Columnar tracks strongly pin vortices, even 
when the tracks are misaligned with the magnetic field by up 
to 40° 06ll . Our imaging work could miss incomplete or mis- 
aligned columns which would be detected in a sample's mag- 
netic behavior. 

For these measurements we use bismuth Hall probes with 
active area 10 x 15 /im^. We place the sample flat atop the 
substrate for the Hall probe. The probe measures total field at 
a spot near the surface of the sample, from which we extract 
the local magnetization of the superconductor We typically 
ramp the applied magnetic field at 15 gauss per minute, and 
take hysteresis loops large enough to allow for full penetra- 
tion. The width of the hysteresis loop, AM, is proportional 
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FIG. 3: Comparison of magnetization and heat capacity for unirradi- 
ated Uo.grTho.osBeis. The vertical lines mark the two transitions. 



to the critical current jc, and we measure both its temperature 
and its field dependence. 

Figure |3l compares hysteresis loop width and heat capac- 
ity for a sample. As T increases, the hysteresis loop width 
decreases, closing at the same superconducting Tc measured 
with heat capacity. For the 3 T thoriated and 7 T pure samples, 
the Tc values in Table HI] are determined from magnetization 
rather than heat capacity data. 

We note that AM{T) in Figure |3 shows no feature at the 
lower transition. By contrast, magnetic relaxation measure- 
ments show a sharp drop in the vortex creep rate at the lower 
transition |6]. UPts, the other HF compound with two transi- 
tions, shows the same behavior: a drop in the relaxation rate 
at the lower transition |6], but no feature in the magnetization 
itself 1 27]. One explanation of the difference between the two 
measurements is that the mechanism of vortex motion changes 
at the lower transition, while the pinning strength itself under- 
goes no sharp change. However, another possibility lies in the 
magnetic field history of the relaxation measurements. Vor- 
tex creep is only suppressed when the sample is exposed to a 
sufficiently large magnetic field, typically 1 kOe, before hav- 
ing the field reduced to zero for the relaxation measurements 
I2al . Our hysteresis loops, with maximum width of 200 Oe, 
may simply be in a different regime where the vortex creep 
shows no signature at the lower transition. 

The magnetization data confirm that columnar defects are 
either absent or far smaller than expected. Increased pinning 
from the columns would appear as wider hysteresis loops, but 
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our loop widths are comparable for unirradiated and irradiated 
samples. In other superconductors columnar defects produce 
large effects from commensurability near integral multiples 
of _B$, but our hysteresis loops have no features whatsoever 
in B. One possibility is that any defects are too much smaller 
than 100 A, the superconducting coherence length iil4ll . to pin 
vortices. 

Our magnetization measurements do show evidence of de- 
fect creation. Figure |3 shows loop width versus normal- 
ized temperature for (U,Th)Bei3 crystals of several irradiation 
doses, taken at zero nominal field. The width follows a power 
law AM (T) = AM (0) (1 - over the entire temperature 
range measured, approximately O.ST^ to T^. For the unirradi- 
ated and the 5 T samples a = 1.8, while for the 7 T specimen 
a is reduced to 1 .4. Although the magnitude agreement be- 
tween the 5 T and unirradiated samples is partly coincidental, 
the 7 T sample has significantly the smallest AM observed. 
Both the magnitude and slope change are probably due to the 
damage. As with the stronger Tc suppression in the 7 T sam- 
ple mentioned earlier, this hints that the 5 T irradiation is a 
threshold for strong damage in our (U,Th)Bei3 crystals. 
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FIG. 4: Hysteresis loop width vs. reduced temperature for three tho- 
riated samples: unirradiated (•), and irradiated with matching fields 
of 5T(A)and7T(D). 



Models for networks of superconducting grains give a = 
1 for superconductor-insulator-superconductor (SIS) tunnel 
junctions |291 and a = 2 for superconductor-normal- 
superconductor (SNS) junctions 1301. An intermediate power 
law for jc{T), with a near 1.5, has been measured in vari- 
ous materials, including the HF material CeCu2Si2 fsT). The 
intermediate values may come from mixtures of SIS and SNS 
junctions in the samples. Similarly, the change in exponent for 
our 7 T sample is consistent with irradiation damage chang- 
ing the characteristics of the grain network. We find the 
same power-law behavior in external magnetic fields up to 
fioH = 2 kG, with the exponent independent of field. 



VI. CONCLUSION 



We have irradiated (U,Th)Bei3 crystals with high-energy 
U and Pb ions at fluences up to _B$ = 7 T. The irradiation 
damage reduces both superconducting transition temperatures 
of Uo.grTho.osBeis by over 15% for the 7 T bombardment, as 
well as lowering jc and reducing its temperature dependence. 
Unlike in previous work with substitutional impurities |23|], 
irradiation causes a similar Tc reduction for the two transitions 
of Uo.97Tho.o3Bei3. 

Despite the measured Tc reduction, transmission electron 
microscope and hysteresis measurements suggest that the de- 
fects are not the amorphous tracks found in HTS. Our calcu- 
lations of the irradiation process predict columnar defect for- 
mation for Se down to 23 keV/nm. Our samples are irradiated 
above the threshold value, at Se ~ 30 keV/nm. The resulting 
defects should have diameter 5 to 6 nm, a defect size com- 
parable to the 10 nm superconducting coherence length. The 
lack of pinning by the defects and their absence from the TEM 
images means either that the actual radius is much smaller or 
that only incomplete columns or point defects form. We favor 
the third option. We rule out the first possibility because of 
the similarity between 208p^56+ 238^^7+ irradiation: if 
the column radius were much smaller for 238^67+ jqjis, then 
no tracks should appear at all for 208p^56+ j^ris. The lack of 
signal in the magnetization indicates that the sample does not 
have even incomplete columns, which would still pin vortices 
effectively. Since the reductions of Tc and jc show that some 
damage is present, we assume we have instead created point 
defects. 

The conditions for creation of columnar defects are im- 
portant to understand, particularly as the use of such defects 
spreads from HTS to other materials. It is unclear whether 
the discrepancy between the experiment and the thermal spike 
calculation comes from the values chosen for the parameters 
or from a deeper problem with the model. The thermal spike 
model has successfully predicted radiation damage both in 
other alloys |9] and in elements, including uranium and beryl- 
lium 1 8]. The simplest explanation would be incorrect pa- 
rameters. Although we use conservative choices, such as as- 
suming resistivity constant rather than linear in temperature, 
our high-temperature extrapolations for p or Kp may be incor- 
rect. Many-electron effects produce a variety of unusual low- 
temperature effects in UBei3 and other heavy fermion materi- 
als, and perhaps could influence higher-temperature behavior 
as well. 

We also show that the lower transition of Uo.gyTho.oaBeis 
is not detectable in jc{T), despite the sharp change in mag- 
netic relaxation rate previously observed. This may arise from 
differences in measurement procedure between our work and 
the relaxation rate studies. Another possibility is that vortices 
in the two superconducting phases of Uo.grTho. 036213 differ 
only in their behavior away from the critical current jc, which 
magnetization measurements would not probe. 
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